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A minimal description of morphological hierarchy in two-dimensional aggregates 

Tamoghna Das,'* T. Lookman,^ and M. M. Bandi*'^ 


A dimensionless parameter A is proposed to describe a hierarchy of morphologies in two-dimensional {2D) aggregates formed due to varying 
competition between short-range attraction and long-range repulsion. Structural transitions from finite non-compact to compact to percolated 
structures are observed in the configurations simulated by molecular dynamics at a constant temperature and density. Configurational ran¬ 
domness across the transition, measured by the two-body excess entropy S 2 , exhibits data collapse with the average potential energy S’ of the 
systems. Independent master curves are presented among S 2 , the reduced second virial coefficient and A, justifying this minimal description. 
This work lays out a coherent basis for the study of 2D aggregate morphologies relevant to diverse nano- and bio-processes. 


1 Introduction 

Aggregates, due to their finite spatial correlations, are struc¬ 
turally intermediate between completely random (liquid) and 
ordered (crystalline) states. Formation of these aggregates 
requires competition^®] between very short range attraction 
and long-range repulsion together with suitable thermody¬ 
namic conditions (low temperature and density). The struc¬ 
ture of aggregates may assume different shapes with varied 
degrees of randomness controlled by the competition® This 
self-assembly of particles is generic and has been observed in 
systems including polymer-coated colloids^®, globular pro¬ 
teins in weakly polar solvents®®, quantum dots and nano- 
particlesESHUl. In addition, the competition can be tuned in 
several ways for a specific system which is advantageous 
from the perspective of industrial applications. For example, 
the range of attraction in a colloidal system can be modified 
by suitable chemical treatment or choice of polymer coating, 
whereas changing the salt concentration or pH balance of sol¬ 
vent controls the effect of repulsion in the same system®!, xhe 
diversity of systems and the flexibility of control have, how¬ 
ever, barred a generic account of aggregate morphologies re¬ 
quired for further pr^ress in developing nano- and bio-metric 
functional materials®!®]. Here we propose a minimal descrip¬ 
tion based on a single dimensionless parameter. 

The descriptor A is a ratio of two effective lengths com¬ 
puted from the competing parts of pair interactions. Tuning 
the interactions, a set of conflgurations is simulated by using 
molecular dynamics. A structural hierarchy within the gener¬ 
ated conflgurations is then characterised by employing stan¬ 
dard statistical diagnostics. Next, the degree of positional ran¬ 
domness governing the hierarchy is quantifled by two-body 
excess entropy S 2 . This enables us to describe all observed 
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morphologies in terms of the average potential energy S and 
A, independently. Collapse of the reduced second virial 
computed for all morphologies, onto a master curve as a func¬ 
tion of A provides experimental access to this description. The 
organisation of the rest of the paper is as follows. In section 
2 we introduce the competing interactions and their control 
parameters including the morphological variations that result 
from tuning the competition with a brief description of simu¬ 
lation details. The spatial correlations for the simulated con¬ 
flgurations and several statistical characterisations of the same 
are presented in section 3. The descriptor A is introduced in 
section 4. In section 5, we quantify the structural randomness 
of observed aggregates and relate it with their respective mean 
energy. We discuss further implications of our flndings in sec¬ 
tion 6 followed by a brief concluding summary. 

2 Competing interactions and variations in 
morphology 

We consider a 2D system of interacting mono-disperse parti¬ 
cles at flxed density p and temperature T. The particles inter¬ 
act pair-wise via a short-range attraction, modelled by gener¬ 
alised Lennard-Jones potential®!, 

0SA=4£[(cr/r)^“-(cT/r)‘*]. (1) 

The length and energy scales of the problem are set by a and 
£ respectively. The range of attraction can be tuned from a 
few times of a to an arbitrary small fraction of it by increas¬ 
ing a. For example, the range is 2.5(7 and 0.2(7 for a = 6 and 
a = 18, respectively. However, the thermodynamics of the sys¬ 
tem does not change for (X > ISEII. The long-range repulsion 
experienced by the particles from their surrounding media is 
considered implicitly by two-body Yukawa potential, 

(/»Lff = (Aa/r)exp(-?-/(*) (2) 

The strength of repulsion A and the screening length ^ is ex¬ 
pressed in units of £ and cr, respectively. The effective centro- 






Fig. 1 (color online) (a) Variation of the effective potential profile (j) (defined in text) with ^ is plotted for A = 4.0 and a = 18. Note that the 
attractive minimum goes from global to local as a function of increasing ^. Inset shows the attractive (jtsA and repulsive (jtiR part of (j) 
separately for ^ =0.5 and same values of a,A. Representative morphologies are shown for ^ = (b) 0.8, (c) 0.7, (d) 0.6, and (e) 0.5, 
respectively. While strong repulsion results in highly anisotropic string-like clusters, the morphology changes continuously to more compact 
shapes as attraction wins over. Further dominance of attraction leads to spanning clusters or gels. Only a part (1/36^^) of the simulation box is 
shown for visual clarity. 


symmetric potential (j) = + (I>lr (Fig (3^) bears an at¬ 

tractive minimum in between a steep hard-core like repulsion 
and a finite positive repulsive barrier decaying exponentially 
to zero at large r. {a, A, } are then three independent param¬ 
eters that tune (j). The variations of (j) is shown in Figjja) as 
a function of ^ for fixed values of A = 4.0 and a = 18. We 
note that the global attractive minimum for =0.5 becomes 
a local minimum and comparable with long-range repulsive 
minimum for ^ =0.8. Varying A for a fixed value of ^ and 
suitably chosen a yield similar variations of (j). The results 
presented next are for a = 18 unless otherwise stated. This 
choice of demonstration is, albeit arbitrary, a reasonable am 
proximation for depletion attraction in globular proteinsE3l31 
and some polymer-coated colloidsE3E3. 

A system of fixed density p =0.4 is used with 56000 parti¬ 
cles in a 376.0(7 x 372.0(7 box with periodic boundary con¬ 
ditions along all directions. Time is measured in units of 
T = for unit mass. Starting from a random configu¬ 

ration of particles equilibrated at high temperature 7^ = 1.0, 
particle trajectories are generated using molecular dynamics 
(implemented by LAMMPSl^S) as the system temperature is 
linearly ramped down to Tf = 0.05 over a period of time 
lO^T. Temperature is measured in £ units and maintained by 
a Langevin thermostat®. The equation of motion for the i-th 
particle with position vector reads as, 

r; = -^V(/»(?-)-vr, + Ci (3) 

considering the force due to interaction (^(r) and frictional 
drag vti from implicit media experienced by the particle. Q 


is a random force with zero mean and Gaussian variance, 
= 2kB{T/v)5ij5{t). Boltzmann constant, 
ks, is chosen to be unity throughout the calculations. Numer¬ 
ical integration of the equation of motion is performed by us¬ 
ing velocity Verlet algorithm with time steps 5t = 10“^ T. We 
specify that our choice of Ty and p is lower than the critical 
values Tc = 0.18 ± 0.01 and Pc = 0.6 ±0.1 for purely attrac¬ 
tive systems®. Aggregation sets in during cooling as soon 
as the system temperature goes below Tc. The linear cool¬ 
ing protocol adapted in this work is very slow compared to 
the typical diffusion timescale (^ t) of the system and is in 
contrast to the quench protocol traditionally followed in ex¬ 
periments. With Tf = 0.05 (standard deviation ^ 10“"^), we 
analyse only the part of the trajectories where the number 
of aggregates does not change over the observation period 
(IO^t). The potential energy, however, continues to decay 
non-exponentially over this observation time. Quantification 
of such non-ergodic behaviour and other anomalous features 
of the microscopic dynamics of the system are described else- 
whereEH in detail. Here, our main focus is on the local struc¬ 
tural properties of the system and its variation with competing 
interactions. Highly anisotropic non-compact clusters of lat¬ 
eral width (7 (Figj^b)) are observed for strong repulsion over 
attraction (large {A,(^}). As a function of increasing attrac¬ 
tion, with decreasing A and/or ^, cluster width increases by 
sidewise aggregation of particles and compact crystalline is¬ 
lands appear in increasing sizes (Fig[2c)-(e)). For even higher 
attraction and near negligible repulsion (small {A,(^}), the 
system consists of one large spanning cluster and several tiny 
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ones. Such states are often referred to as gel^^. 

3 Structural characterisation 

3.1 Radial and angular correlations 
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Fig. 2 (color online) (a) Radial distribution g{r) of particles for the 
same configurations shown in Figj^h) and (e). While compact 
structures =0.5) show symmetric splitting of second peak 
characteristic of crystalline order in the system, this feature becomes 
considerably less and anisotropic for non-compact aggregates 

=0.8). (b) The distribution of three-body angle P( 03 ) is 
computed for the same systems as above. In accord with previous 
observations, P( 03 ) shows three preferred peaks at 7r/3,27r/3 and % 
for compact crystalline arrangement. Illustrations over each peak are 
representative of particle arrangements for the three different cases. 
Last two peaks (03 =27r/3, tt) become less prominent for 
non-compact structures where particles arrange in a variety of 
angles to each other, for example, see Figl^a). 

To characterise the diverse morphologies observed, we now 
focus on the microscopic features of particle arrangements in 
simulated configurations. The radial distribution function® 

g{r) = \/p(^5{r-ri)) (4) 

computes the probability of finding an i-th particle from an 
arbitrary central one as a function of the separation between 
them. Whereas the first peak in g(r) (Fig|^a)) accounts for the 
nearest neighbours, splitting of the second peak indicates two 
preferential next nearest neighbour distances, typical of crys¬ 
talline (compact) arrangements (Figjje)). For non-compact 
aggregates, the peak heights reduce and the second peak be¬ 
comes anisotropic. The long shoulder between the preferential 


positions of first and second peaks accounts for the possible 
non-crystalline arrangement of particles observed in Figj^b). 
g{r) for other intermediate structures (Fig[^c)-(d)) naturally 
falls within these two limits with varied peak heights and 
anisotropy accordingly. Angular arrangement of particles, be¬ 
ing inaccessible by g{r), is studied using three-body angles 03 . 
All triads of particles are considered where each particle is the 
nearest neighbour of at least one other particle and the dis¬ 
tribution of three-body angles PiO^) (Fig|^b)) is computed. 
Three preferential values, ;r/3,2;r/3 and 7i, assumed by O 3 are 
indicative of hexagonal symmetry expected in 2D. Whereas 
compact aggregates agree with this, non-compact ones nat¬ 
urally deviate as the hexagonal symmetry breaks down and 
more angular arrangements become possible. The radial and 
angular features of all percolated configurations are similar to 
the compact ones, though enhanced as expected. 

3.2 Shape and size statistics 



Fig. 3 (color online) (a) Fractal dimension df of aggregates is 
computed from the slope of their size 5 plotted against radius of 
gyration Rg. Here, we plot the values computed for the 
representative configurations shown in Fig[2b)-(e). Long- and 
short-dashed lines have slopes 2 and 1 respectively, (b) Size 
distribution P{s) of aggregates for different values of ^ at A = 4.0. 
Power-law behaviour of P{s) is indicative of percolation at ^ =0.4 
in contrast with the exponential behaviour for other values. A 
dashed line of slope —2 is shown for comparison. 

Visually evident structural transition of finite-size aggre¬ 
gates is further characterised by their fractal dimension df. 
Exploiting the relation between the size 5- of a cluster and its 
radius of gyration Rg, 

s ~ (5) 

R] = ^/'^sY,{ri-rjf 

J/ is obtained for different conformations. For the non¬ 
compact (almost linear) and compact (locally crystalline) ag- 
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gregates, J/ is bounded between 1 and 2 corresponding to lin¬ 
ear and planar shapes respectively (Fig|^a)). As the transi¬ 
tion occurs seamlessly with varying competition, it is difficult 
to define a sharp boundary of structural transition in terms 
of df. We consider all configurations with Jy > 1.6 as com¬ 
pact and the rest as non-compact. All these finite-size ag¬ 
gregates appear in exponentially distributed sizes within the 
system. These clusters aggregate further as attraction domi¬ 
nates over repulsion and span the system as a single perco¬ 
lating cluster. The size distribution P[s) for the ensemble of 
clusters thus charges from exponential to algebraic and be¬ 
comes scale-free® 


P{s)^s ^exp(^/^o) (6) 

where v is the Fisher exponent and exponential part is due to 
the finite-size effect (Figj^b)). Fitting the cluster size distri¬ 
bution with the above mentioned form, we numerically esti¬ 
mate the Fisher exponent V = 1.98 ± 0.04. This value is very 
close to the random percolation (RP) model Fisher exponent, 
187/91, in 2D. Within RP model, V < 2 is spurious as it im¬ 
plies unphysical divergence of mean cluster size. However, 
v = 1.91±0.06 is reported experimentallyEl and can be ex¬ 
plained through a simple extension®! of RP model. Identi¬ 
fication and characterisation of the exact nature of this geo¬ 
metric transition require further in-depth investigation and we 
leave that for future correspondence. We mention that the sys¬ 
tem density (p = 0.4) is lower than the typical critical value. 
Pc = 0.5, for the RP model. The exponential cut-off is only 
very weakly dependent on repulsion controlled by {A, and 
increases with density. Keeping the interaction intact, similar 
percolation has also been observed as a function of increas¬ 
ing density. The non-compact aggregates percolate at p =0.5, 
consistent with typical random percolation. 


4 Single descriptor for aggregates’ morphology 
phase diagram 

We now compute an effective hard-core diametei®®!, 

Oi= [ [1-exp(-j30/,)]Jr (7) 

Jo 


set by (j)sA where 0/^ = 4e(cr/r)^^. cJi is thus a function of £ 
and a. Two control parameters of repulsion, namely, A and ^ 
can be similarly encoded into another length scale. 


02 = 


poo 

/ [^-ew{-P(l>LR)]dr 

Jo 


( 8 ) 


where j3 = 1/ {ksT) with Boltzmann constant, = 1. As the 
definitions are not specific to the functional form of the poten¬ 
tial, they can be computed for other forms®!S2 

of competing 


interactions with proper care. For example, it is enough to 
compute only G 2 for 2-Yukawc^ov Coulomb-YukawJ^ mod¬ 
els as the extent of hard-disk interaction is fixed in such mod¬ 
els to match the real physical situations at hand. The effective 
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Fig. 4 (color online) (a) The dimensionless parameter A (defined in 
text) is plotted against ^ for the same set of A’s with two different 
a= 18 (open symbol) and 12 (filled symbols). Please note that one 
specific value of A can be obtained through several combinations of 
{a,A, (^}. (b) Morphology phase diagram: boundaries for different 
structures drawn on the map of A in A-(^ plane for a = 18. 
Non-compact to compact transition determined by the fractal 
dimension of respective aggregates is marked by the dotted line; the 
dashed line denotes percolation transition determined from the size 
statistics of clusters. 


lengths, a I and (72, project a useful illustration of the system : 
consider a system of hard particles with effective diameter a \, 
each having a concentric soft shell of influence circle with di¬ 
ameter c 72. Inclusion of particles within G 2 is possible which in 
turn changes the perimeter of infiuence zone of the combined 
particle system in a non-trivial way thus determining the ag¬ 
gregate morphology. The ratio, A = ( 72 / (7i, then encodes all 
three control parameters, a, A and ^ and several combina¬ 
tions of these parameters result in the same A. Fig^a) points 
out that small a-large A is equivalent to large a-small A for a 
range of (A for a = 18, A=3.0 almost coincides with a = 12, 
A = 5.0). Larger a(> 18) will populate the upper left half of 
Figga), whereas smaller a(< 12) will belong to the lower 
right half of the same for the specified rang e of A and ^. The 
hard disk limit with strict exclusion region®MI recovered 
for A = 1. Mapping of {a, A, } space to a single parameter A 
indicates that the morphologies controlled by the former pa¬ 
rameters can possibly be represented by the latter. Now, we 
draw the boundaries of different conformations on a A-map of 
{A, parameter space for fixed a = 18 (Fig|^b)). This will 
serve as the morphology phase diagram of the present model 
system. Following the appearance of power-law size distribu¬ 
tion of clusters, we find that percolated structures or gels are 
always formed for parameters giving A< 1.5. For larger A, the 
aggregates are finite sized and change from compact to non¬ 
compact shape around A ^2.0 which can be verified from the 
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fractal dimension df of respective aggregates.. Taken together 
with our previous observation, we suggest that lowering a will 
lead to compact and gel structures whereas non-compact ones 
are favoured by large a for this system. 

5 Structural randomness and mean energy 



Fig. 5 (color online) Two-body excess entropy S 2 provides a 
quantification of positional information. Collapse of ^2 with A 
offers a unified description of aggregate morphologies independent 
of their control parameters {a,A,^}. Comparing with morphology 
phase diagram (FigQb)) for a = 18, we can now identify 
non-compact (A>2.0), compact (1.5 < A<2.0) aggregates and gels 
(A< 1.5). The dashed and dotted lines are phase boundaries used in 
Figj^b). As shown, this mapping also holds for a = 12. 


How are the observed morphologies related to their ener¬ 
gies? Answering this question requires further quantification 
of particle arrangements. Positional information of particles 
can be appropriately expressed as multi-particle expansion of 
the total excess entropyl^ (relative to an ideal gas of same 
density) truncated at the two-body term, 

52 =-p/2 j (/r{ 5 (r)lng(r)-[g(r)-l]}. (9) 

This ensemble invariant measure® of disorder is 0 for maxi¬ 
mal disorder (ideal gas) and goes to — oo for the ordered (crys¬ 
talline) state. When plotted against A (Fig0, ^2 shows three 
distinct ranges corresponding to three different morphologies 
charted out qualitatively in the previous morphology phase 
diagram (Fig^b)). The non-compact structures for large A 
posses higher degrees of randomness close to the random con¬ 
figurations. We note that these structures appear under strong 
influence of long-range repulsion. One particle can, however, 
come within the large repulsive influence zone (02 > 2cri) of 


another particle and is unable to escape due to short-range at¬ 
traction. This situation evidently favours a collective unidi¬ 
rectional arrangement and opposes the centro-symmetric fea¬ 
ture of the effective pair-wise interaction. As the repulsion 
decreases, particles can not be accommodated within each 
other’s repulsive influence zone and the local hexagonal sym¬ 
metry favoured in 2D is recovered. S 2 reflects this feature by 
showing a sharp drop for A <2. As G2 and G\ becomes compa¬ 
rable (A< 1.5) under strong influence of attraction, percolated 
conformations appear and ^2 assumes even larger negative val¬ 
ues. Collapse of ^ 2 , computed for conflgurations with a = 12 
(Fig|^ on the same curve as a = 18, establishes the robustness 
of this mapping against the attraction range. The morphologi¬ 
cal hierarchy of 2D aggregates can then be described in terms 
of A alone, and independent of the individual control parame¬ 
ters 



Fig. 6 (color online) (a) Average energy i versus A is linear; the 
slopes increase with A. (b) Plotted against average energy S, S 2 
clearly shows three distinct ranges. These ranges with decreasing 
energy account for three candidates of morphological hierarchy, 
non-compact, compact and percolated structures, respectively. 

Since the thermal fluctuations are very small at low temper¬ 
atures, potential energy plays the dominant role in determining 
local structures. The average potential energy # is deflned as: 

# = (1 / A) Ei where Et is the potential energy of i-th particle 
due to all other particles present within a cut-off radius of lOcr. 

# follows a linear relation with A for flxed A and the slope 
steepens with increasing A (Fig|^a)). In other words, as the 
soft shell of repulsion increases, it is possible to And a set of 
degenerate A’s with same energy and vice versa. As measured 
by S 2 , the morphological hierarchies can now be identifled by 
their distinct ranges of average energy # (Figgb)). Highly 
disordered non-compact structures (S 2 ^ 0) are spread over a 
range of #. In contrast, compact structures occurring with a 
small variance of # show large variation of S 2 . Gels have even 
lower values and smaller spread of both ^2 and #. These ob¬ 
servations certainly call for an in depth understanding which 
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we leave for future investigation. 
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6 Discussion 

We note that the knowledge of interactions is necessary for 
this description. By relating B 2 , the second virial coefficient, 
with A, we relax this requirement and connect the description 
with experiments. Experimental estimation of B 2 does not re¬ 
quire any prior knowledge of specific potential and/or confor¬ 
mation for a real system. For a simulation study, however, B 2 
is computed as followsE^, 

B 2 = -{1/2T) Jdr<l>\r)g{r) (10) 

using both potential and radial distribution. The prime over 
(j) denotes the first order spatial derivative of the same. The 
reduced second virial coefficient, B\ = 2 B 2 /, showing 
a spread as a function of A and (^, follows a master curve 
when plotted against A (Fig|^ for both a = 18 and 12. We 
note that most experimental efforts specify the overall nature 
of interactions but do not probe the functional form that cor¬ 
responds to a given measurement. While the morphologies 
have been described in several ways, g(r) is, in general, not 
readily available to calculate B\. For example, the existing B\ 
data for globular proteins is inadequate to specify the interac¬ 
tions and/or structures completely. The B\-K master relation 
provides access to the microscopic length scales from thermo¬ 
dynamic measurements. This collapse has further conceptual 
implications. B\ has already been successfully used^SIto find 
thermodynamic correspondence among wid e rang e of systems 
from the van der Waals limit (simple liquids)E313Ito the Baxter 
limit (sticky hard spheres)ISl. This study shows that B\ is ca¬ 
pable of capturing the morphological correspondence as well 
and providing a feasible basis for the extension of correspond¬ 
ing states principl^^^io aggregates. _ 

A formal order parameter approach, mean fieldE^SI 
beyondlSl is often challenging for pattern forming systems 
showing wide variations in local density. Quantification of 
positional randomness by ^2 provides an alternative descrip¬ 
tion capable of relating the local geometry of conformations 
to a thermodynamic quantity like #. Efforts have been made 
to describe the aggregates employing the integral theory ap- 
proach^ZIfor simple liquids. However, such an approach con¬ 
sidering the effect of repulsion as a perturbation to attraction 
is inadequate for the systems presented here where both of the 
competing interactions are comparable in strength. We have 
been able to encode the details of competing interactions into 
two distinct length scales using standard liquid states theory. 
The procedure is valid over a wide range of control parameters 
and should also be readily applicable to other model systems 
as pointed out earlier. Further, these length scales offer an 



12 


8 


4 


0 


Fig. 7 (color online) Second viral coefficient B^ computed for the 
full {a,A, parameter space collapses onto a master curve when 
plotted against A. Same symbols as in Figj^have been used. Inset 
Please note the spread in when plotted as a function of ^ for 
three different values of A for a = 18. 

intuitive grip over the physical situation useful for further the¬ 
oretical developments. The ratio of these length scales, A can 
successfully identify all three major candidates, non-compact, 
compact and percolated structures, forming the morphological 
hierarchy of aggregates in 2D. We mention that availability of 
more degrees of freedom in three dimensions may potentially 
affect the scenario presented here. However, the framework 
should be valid and we leave the detailed validation for fu¬ 
ture investigation. Importantly, A is experimentally accessible 
through the master relation between ^nd A. Given the 5'2-A 
map, the expected morphologies can then be easily predicted. 
This can be potentially important for the design of new ma¬ 
terials with desired functionalities. In conclusion, we have 
presented a minimal, yet robust description of 2D aggregate 
morphologies potentially relevant for the self-organisation of 
diverse nano- and bio-processes. The present work provides 
the missing information connecting interactions, microscopic 
structures and thermodynamics. This should motivate further 
studies on complex fluids probing their structural signatures 
with higher resolution. 
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